Study of Hund's rule coupling in models of magnetic impurities and quantum dots 
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Studies of the effects of the Hund's rule coupling Jh in multiple orbit impurities or quantum 
dots using different models have led to quite different predictions for the Kondo temperature 7k 
as a function of Jh. We show that the differences depend on whether or not the models conserve 
orbital angular momentum about the impurity site. Using numerical renormalization group (NRG) 
calculations, we deduce the renormalized parameters for the Fermi liquid regime, and show that, 
despite the differences between the models, the low energy fixed point in the strong correlation 
regime is universal with a single energy scale 7k, and just two renormalized interaction parameters, 
a renormalized single orbital term, U = 47k, and renormalized Hund's rule term, Jh = 87k/3. 



I. INTRODUCTION 

The role of the Hund's rule coupling Jh has been 
investigated for magnetic impurities and quantum dots 
by a number of authors starting from rather different 
models— 5 . This has led to different predictions for the 
behavior of the Kondo temperature in these systems as 
a function of Jh- As a result it is not clear whether or 
not the models may also differ in their predictions for 
their low energy behavior. We investigate this question 
by applying a combination of the numerical renormaliza- 
tion group (NRG) and renormalized perturbation theory 
(RPT) to some of these different models. From an anal- 
ysis of single particle excitations about the NRG low en- 
ergy fixed point we can obtain an accurate estimate of 
the Kondo temperature in all cases. We can also calcu- 
late the renormalized interaction parameters at this fixed 
point and compare them for the different models. This 
information can be used to deduce the low temperature 
specific heat coefficient, spin and charge susceptibilities 
and low energy dynamical behavior for the models in all 
parameter regimes. 

A model of a magnetic impurity with n-fold degenerate 
orbitals interaction with a bath of conduction electrons 
will in general require quite a large number of parame- 
ters to describe the interactions between the electrons in 
these orbitals for various filling factors, and their interac- 
tion with the conduction electrons of the host metal. To 
understand the basic physics of these models a number 
of different models, using a restricted set of parameters, 
have been used. One of these models was introduced 
by Yoshimori (model 1) where the interactions between 
electrons in the impurity orbitals were specified by just 
two parameters-, a local Coulomb interaction U and an 
inter-orbital exchange interaction Jh. The Hamiltonian 
of the model takes the form, 
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gular momentum quantum number I, and z-component 
to = —I,— I + 1, • ■ • ,1, and spin component a =t,i- 
The impurity level in a magnetic field H we take as 
Edma = £rf - MBfiJ - fi B mH - [i, where er = 1 (f) and 
(7 = —1 (J.) and /i is the chemical potential, and hb 
the Bohr magneton. The creation and annihilation op- 
erators ct , c kma are for partial wave conduction elec- 
trons with energy e^ma- The hybridization matrix el- 
ement for impurity levels with the conduction electron 
states is Vk- We denote the hybridization width factor 
by A mCT (e) = 7r^] fc \Vk\ 2 S{e-ekma), which we can take to 
be a constant A in the wide flat band limit. The remain- 
ing part of the Hamiltonian, rid describes the interaction 
between the electrons in the impurity state, 
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where d\ na , d ma , are creation and annihilation opera- 
tors for an electron in an impurity state with total an- 



This model can be used to describe transition metal im- 
purities, such as Fe or Mn, in a metallic host in the ab- 
sence of spin orbit or crystal field splittings. The model 
can be interpreted more generally with a = to + I + 1 
as a channel index taking values a — 1, 2, • • • , n where 
n = 21 + 1 is the number of channels. The Hund's rule 
term tends to align the electrons on the impurity site 
such that for large U and large Jh the impurity state 
will correspond to a spin S — n/2. 

A different n-fold model has been studied recently by 
Nevidomskyy and Coleman 2 (model 2) which also in- 
cludes a Hund's rule interaction. This model is char- 
acterized by just two interaction terms, an inter-orbital 
Hund's rule exchange interaction Jh and an Kondo form 
of exchange interaction Jk between the impurity elec- 
trons and the conduction electrons, so that the Hamilto- 
nian takes the form, 

H mo dc\2 — J2k,aa e ka<jc\ aiy C karJ — Jii(J2a=l Sq ) 2 

+ J K E«=l Efc,fc' Sa-Cfcaa^'CfcW ( 4 ) 

There are significant differences between models 1 and 2, 
which can be seen more clearly if we restrict the discus- 
sion to the case n — 2. when model 1 can be written in 
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a more usual form with the exchange interaction written 
in terms of spin operators, 



Hd. 



model 3 
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(5) 

with U12 = U— 3Jh/2. We generalize the model by taking 
U\2 as an independent parameter, and refer to the model 
in this generalized form as model 3. It has been used in 
this form to describe the interactions in certain double 
quantum dots^. 

We consider, first of all, model 3 in the limit Jh = 
with U\2 = 0. In this case there is no interaction or 
hybridization between the two channels and the model 
is equivalent to two independent Anderson models. In 
the strong coupling regime U / ir A 3> 1 with particle-hole 
symmetry it can be converted into model 2 ( Jh = 0) via a 
Schrieffer- Wolff transformation with Jk = 4 V 2 /U. If the 
first order correction term in Jh is included in this limit, 
then model 3 can be seen to be essentially equivalent to 
model 2 in Eq. (@J for n = 2. 

However, for Jh = 0, and U12 = U in Eq. ([5]), corre- 
sponding to the Yoshimori model (model 1), there is an 
interaction between the two channels. In this case the 
orbital index I can be combined with the spin index into 
a single index (to, a) running over 4 values and the model 
can be shown to have SU(4) symmetry, corresponding to 
conservation of both spin and orbital angular momen- 
tum. In the particle- hole localized limit with U / ir A 3> 1 , 
it can be transformed via a Schrieffer- Wolff transforma- 
tion (which takes into account virtual excitations from 
the 2-electron, six-fold degenerate, ground state of the 
impurity to 1-electron and 3-electron excited states) into 
an SU(4) Kondo model with J K = W 2 /U with a six di- 
mensional representation for the local SU(4) operators 3 . 
Hence models 1 and 2 are quite different in the Kondo 
regime for Jh = 0. For general n, in the Kondo regime 
when Jh = 0, model 1 reduces to an SU(2n) Kondo 
model and model 2 to n independent SU(2) Kondo mod- 
els. This reflects the important difference that model 1 is 
invariant under orbital as well as spin rotation, whereas 
model 2, and model 3 when U\ 2 ^ U — 3Jh/2, are in- 
variant under spin rotation only. This can be verified ex- 
plicitly by considering the commutator of the generator, 



+ c 
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2(T ), with the Hamiltonian 



of model 3. This difference could have consequences for 
the low energy behavior of these models which we con- 
sider in the next section. 



II. LOW ENERGY FERMI LIQUID REGIME 

The conservation of both angular momentum and spin 
in model 1 permitted Yoshimori^ to derive two Ward 
identities for interaction vertex part at zero frequency. 
This interaction vertex can be interpretated in terms 
of a renormalized interaction between the quasiparti- 
cles of a Fermi liquid 3 described by the Hamiltonian, 



H = Ho + Hd, where 

m,(T ka 

) (6) 



and 



n d = 



(U - J H ) 
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This effective model is of the same form as the original 
model defined in Eqs. © and Q with the difference that 
the interaction terms have to be normal ordered. The 
brackets : O : indicate the normal ordering of the opera- 
tor O with respect to the ground state of the interacting 
system, which plays the role of the vacuum. This term 
only comes into play when more than one quasiparticle 
is created from the vacuum. 

From the Ward identities 3 -^ exact expressions can be 
derived for the impurity charge, spin and orbital suscep- 
tibilities, Xc Xs, an d Xorb, at T = 0, 

Xc = 2np(°) (0) (1 - ((2n - 1)U - 3(n- I) J H )/5 (0) (0)). (8) 



Xs - 2n M |p (0) (0)(l + {U+{n- I) J H )p (0) (0)), (9) 



X^= ( " a -y (0) (l + (^-3Jk)^(0)), 

(10) 

where p(°\uj) is the free quasiparticle density of states 
per single spin and channel, 



/5<°>( W ) = 
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The impurity contribution to the specific heat coefficient 
7 is also given exactly by 7 = 2n7r 2 p(°)(0)/3 and the 
Wilson ratio i?w = 7r2 Xs/3/^ B 7 is given by 



R w = 1 + (U + J H )/5 (0) (0). 



(12) 



In the localized limit U / ttA 3> 1 with particle-hole 
symmetry (e<j = 0), the charge susceptibility is sup- 
pressed so we can equate Xc to zero. Similarly when 
Jh is large the orbital susceptibility will be suppressed 
so that Xorb can also be equated to zero. Then from Eqs. 
© and (jTUJ) we obtain for the spin susceptibility, 



(g^ B ) 2 S(S+l) 
3T K 



(13) 
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where S = n/2 and g = 2, and the relation between the 
renormalized parameters, 

ttA = U = | J H = 4T K . (14) 

As angular momentum is not conserved for model 2, 
and model 3 if U12 ^ U — 3Jh/2, the question arises as 
to whether the low energy fixed point of this model can 
be described by a quasiparticle Hamiltonian with renor- 
malized parameters, similar to that given in Eqs. ([5]) and 
([7]), and if so, does the relation given in Eq. (fT4"|) hold 
when Jh is large? 

To examine this question we extend our earlier NRG 
calculations- of the renormalized parameters for model 
1 to models 2 and 3 for the case n = 2. We can test 
the hypothesis that the low temperature behaviour of all 
three models can be described by a quasiparticle Hamil- 
tonian in terms a set of renormalized parameters, A, U, 
U12 and Jh- We restrict our calculations to the particle- 
hole symmetric case and take €d = — U/2 — U12 in the 
one-electron part of the Hamiltonian given in Eq. (J2J. 
In the NRG calculations we used for the discretization 
parameter A = 6 and half-bandwidth D = 1, and ap- 
proximately 4000 states were retained at each iteration. 
The value 7rA = 0.01 was used for models 1 and 3. 



III. CALCULATION OF RENORMALIZED 
PARAMETERS AND T — SUSCEPTIBILITIES 

The renormalized parameters that describe the quasi- 
particles and their interactions can be deduced from an 
analysis of the low energy NRG fixed point. The param- 
eters id and A can be deduced by fitting the lowest single 
particle and hole excitations from the NRG ground state 
to those of a non-interacting Anderson model. The inter- 
action parameters, U, Ui 2 and Jh, can be deduced from 
the lowest two-particle excitations, in the same channel 
for U, and different channels for U12 and Jh- We have 
analysed the fixed point in the same way for all the mod- 
els, including model 2. Details of these calculations were 
given in earlier worker— so we shall just quote the re- 
sults of such an analysis here. 

We use the relation, nA = 4Tk to deduce the Kondo 
temperature in the particle-hole symmetric case for all 
the models. The results for 2k(Jh)/Tk(0) as a func- 
tion of Jh/4Tk are shown in Fig. [1] for particular pa- 
rameter sets for the three models. For model 1 with 
U/ttA = 5, the value of Xk(Jh) falls of rapidly with 
increase of Jh- Results for this case were presented 
in earlier work^, and shown to decrease exponentially 
with a fit Tk/ttA = 0.0854exp(-1.497r 2 J h /ttA) with 
7rA = 0.01. This exponential dependence is to be 
expected from a Schrieffer- Wolff transformation in the 
regime when U/ttA 1 and Jh is large, when the 
ground state of the isolated impurity is a two electron 
spin triplet state. When virtual excitations to local states 
with one more or one less electron are taken into this 



leads to a spin 1 Kondo model with an effective antifer- 
romagnetic coupling J£ = W 2 / {U + 3Jh), giving expo- 
nential factor in Jh in the resulting Kondo temperature 
T K ~ De~ 1 ' 2J x pc , where p c is the density of states of the 
conduction electrons at the Fermi level. 

The results for model 2 are for the case with Jk = 0.3, 
and are seen to decrease much more slowly with increase 
of Jh and appear to level off for larger values of Jh- 
However, the initial values of Tk(0) are rather differ- 
ent for the two models. As model 2 in this limit re- 
duces to n-independent SU(2) Kondo models the Kondo 
temperature is given by Tk(0) ~ De~ 1 / 2jKPa . In this 
limit, on the other hand, model 1 corresponds to an 
SU(2n) Kondo model, which has a Kondo temperature 
T K ~ £) e -i/2«J K pc ; so for the same C0U pii n g t k {0) is 

likely to be significantly larger than that for the SU(2) 
model. 

The results for model 3 with U/ttA — 5, U12 — 0, also 
shown in Fig. [TJ are very similar to those of model 2. This 
is not surprising because because, as noted earlier, when 
U\2 = and Jh = it can be transformed into model 2 
(with Jh = 0) via a Schrieffer- Wolff transformation when 
U /ir A is large, and adding the first order correction term 
in Jh gives essentially model 2. 

Another parameter set for model 3 with U/irA = 3, 
U12/1TA = 2.9, is shown in Fig. [TJ In this case there 
is still a marked fall off of Tk with Jh but intermediate 
between the other cases shown. 




y4T K (0) 

FIG. 1. (Color online) A plot of T k (Jh)/T k (0) versus 
Jh/4Tk(0) for model 1 with [7/ttA = 5, model 2 with 
Jk = 0.3, model 3 with (7i2/7rA = 2.9, U/nA = 3 and also 
with Uw/irA = 0, U/irA = 5. 



In Fig. [2] we examine the dependence of Tk on Jh 
for model 2 in more detail by plotting ln(Tk(JH)/Tk(0)) 
versus ln( Jh/4Tk(0)) for a range of values for the Kondo 
coupling, J K = 0.1,0.12,0.15,0.2,0.25,0.3. All the re- 
sults fall on a universal curve provided Jk < Jh- At 
J K ~ Jh the curves begin to deviate and develop a 
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plateau for Jh > Jk for the larger values of Jk, Jk = 
0.15,0.2,0.25,0.3. Plateaus also develop for the smaller 
values of Jk for Jh/^k(0) beyond the range shown in 
Fig. [21 We can compare these results with the scal- 
ing relation derived to one loop order for this model by 
Nevidomskyy and Coleman 2 (Eq. (11) in their paper 
with n = 2) for intermediate values of the Hund's rule 
coupling Jh, which predicts Jk(Jh) °c ?k(0)/Jh- This 
would correspond to a linear behavior in the log-log plot 
in Fig. [5] with a slope —1. The dashed line Fig. [5] cor- 
responding to Tk(Jh) = 4T^(0)/e 2 Jn does fit the curves 
over a substantial range, particularly for the smaller val- 
ues of Jk- However, a linear fit to the range Jh < Jk 
for larger values of Jk would give a larger slope with 
Jk(Jh) oc (Jh) _x , with x varying 1 < x < 1.25 depend- 
ing on the fitting range chosen. The linear fit with x ~ 1 
begins at values of Jh ~ 4Tk(0). 
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FIG. 2. (Color online) The graph shows a plot of 
1ii(Tk(Jh)/Tk(0)) versus m(J H /4T K (0)) for model 2 for val- 
ues of J K = 0.1, 0.12, 0.15, 0.2, 0.25, 0.3. The dashed line cor- 
responds to Tk(Jh) = 4T K (0)/e 2 J H , on using the relation 
4Tk(Jh) = 7tA(Jh). The deviation from the approximate lin- 
ear dependence to a plateau region for Jk = 0.15, 0.2, 0.25, 0.3 
occurs when Jh ~ Jk in all four cases. 



When Jh Jk and the two impurity spins are locked 
together, model 2 should also be equivalent to an effec- 
tive Kondo model with the electrons in the two channels 
coupled to a spin S = 1. For the general n-channel model 
this Hamiltonian would take the form, 



H = Jt 



m=l k,k' 



(15) 



with an effective coupling J K of a spin S = n/2 to the 
conduction electrons. In the limit Jh —> oo the result, 
Jk = Jk/ti, derived by Schrieffer— for model 2 should be 
valid. The Kondo temperature of the model with n = 2S 



would then be given by T K ~ De 1 / 2J kP<> } neglecting any 
dependence of the prefactors on J K . The corresponding 
result for the model with J H = 0, T K (0) ~ De~ l / 2J *p c , 
again neglecting any dependence of the prefactors on Jk, 
would then imply T K oc (T K (0)) n . We can test this result 
for n = 2 by modifying the plot given in Fig. [5] and plot 
instead ln(T K (J H )/T|(0)) versus ln( J H /4T K (0)). The re- 
sults are shown in Fig. [3]for Jk = 0.3, 0.25, 0.2, It can be 
seen for the cases shown that, in the regime Jh > Jk> the 
ratios of Tk(Jh)/?k(0) become independent of the value 
of J K which implies that, for J H > Jk, T k oc (T K (0)) 2 . 
This gives a clear criterion Jh > Jk for the Schrieffer 
relation^, J£ = Jk/^, to hold. 
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FIG. 3. (Color online) A plot of ln(T K ( J H )/T K (0)) versus 
ln(J H /4T K (0)) for model 2 for values of J K = 0.3,0.25,0.2. 
The results indicate that the ratio Tk(Jh)/Tk(0) becomes 
independent of Jk in the regime Jh > Jk- 



In Fig. [4] we plot the ratio of the renormalized param- 
eters, J#(Jh)/7tA(Jh) versus J^/4Jk(0) for the differ- 
ent models for the parameter sets given in Fig. [TJ We 
can see that all the curves asyptotically approach the 
value predicted in Eq. ([11]), Jh(Jh)/VA(J h ) ->• 2/3 on 
increasing Jh- Though they approach this value at dif- 
ferent rates they are all very close to the limiting value 
for Jh > 57k (0). This clearly shows that the crossover 
to T = susceptibility corresponding to a spin 1 Kondo 
model, as given in Eq. (|13p . on increasing Jh occurs on 
a scale Tk(0). The regime 5Tk < Jh < Jk f° r model 
2 corresponds to the scaling regime Tk(Jh) oc 1/Jh of 
Nevidomskyy and Coleman^ where the low temperature 
behaviour can be described by a spin 1 Kondo model be- 
fore the crossover to the regime Jh > Jk where Jk(Jh) 
becomes independent of Jh, and the Schrieffer result 1 
J K = J K /2, T* kT 2 (0), holds. 

The remaining renormalized parameter ratios, 
[//7tA(Jh) and U^/ttA^Jh), that characterize the low 
energy fixed point, are shown in Fig. [5] for the different 
models as a function of Jh/4Tk(0). For model 1 with 
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U/ttA = 5 the charge fluctuations are suppressed Xc ~ 0, 
so from Eq. ([5]) when Jh = (Jh = 0) for n = 2 we 
get JJ/ttA = 1/3. As Ux 2 = U - 3J H /2 for this model 
when Jh = 0, we also predict the result U^/irA = 1/3. 
These are confirmed in the results given in Fig. [S] When 
Jh is increased and the orbital fluctuations are also 
suppressed we get the results for this model given in Eq. 

which imply that in this limit Uyi — > 0, and the 
NRG calculations confirm these results. 

The NRG results for model 2 give U/ttA = 1 and 
U12 = for all values of Jh- We also obtain the same 
results for model 3 (not shown) for the case U\2 = 0, 
U/ttA — 5. This is not surprising because, as noted 
earlier, when [/ 12 — > and U^/itA 1, model 3 and 
model 2 are very similar. 

Finally in Fig. \5\ results are shown for model 3 for 
the case t/^/vrA = 2.9 and U/nA = 3. To test the 
predictions in this case, the expression for \c given in 
Eq. ([5]) must be generalized as it assumes the relation. 
U12 = U — 3Jh/2. As both charge and spin are conserved 
in model 3, from the corresponding Ward indentities, ex- 
act results for the charge and spin susceptibilities can be 
derived. For n = 2 the generalized expression for the 
charge susceptibility takes the form, 

Xc =¥ 0) (0)(l-U-2U 12 )pW(0)). (16) 

As the charge fluctations are suppressed for model 3 with 
Ui 2 /nA = 2.9 and U/ttA = 3, Eq. (H]) predicts the 
result U/ttA + 2J7 12 /7rA = 1, which is satisfied in the 
results shown in Fig. [3] for all values of Jr. Though 
both U /it A and V \ 2 j 'tt A are non-zero for Jh = 0, as Jh 
is increased U /w A — > 1, U^/ttA — > in line with the 
result in Eq. (fT4|) . 




J H /4T K (0) 

FIG. 4. (Color online) A plot of jn/nA(Jn) versus J H /4T K (0) 
for the models and parameter sets given in Fig. [T] 
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FIG. 5. (Color online) A plot of U/ttA{J u ) and U 12 /irA{J H ) 
versus Jh/4Tk(0) for the models and first three data sets 
shown in Fig. [1] The corresponding results for the last data 
set given in Fig. [T]is U/irA(J H ) = 1 and U 12 /-kA(J h ) = for 
all values of Jh in the range shown. 



The expressions for spin susceptibility \s and Wilson 
ratio i?w for model 3 with C/ 12 ^ U — 3Jh/2 remain 
unchanged from those given in Eqs. (j9]) and (|T2|) . To 
calculate the impurity spin susceptibility for the models 
and parameters sets used in the other plots we substi- 
tute the renormalized parameters into Eq. ([9]) and give 
the results in Fig. [51 The exponential decrease of the 
Kondo temperature with increase of Jh gives a dramatic 
rise in \s for model 1 compared with the corresponding 
results for models 2 and 3. However, the change in the 
corresponding Wilson ratios on increasing Jh from zero 
is much less dramatic. For model 1 with U /irA = 5, i?w 
increases from 4/3 for Jh = to 8/3 for large Jh, and 
for model 2 from 2 to 8/3 over the same range. 



Yoshimor: 6 derived an exact result for the low temper- 
ature impurity contribution to the resistivity for model 1 
in the particle-hole symmetric case and H — 0. In terms 
of the renormalized parameters, the result is 

^) = i?o(l-^f^(^) 2 + 0(T^, (17) 

where / is given by 

/ = 2((2n - l)U 2 - 6(n - l)J n (U - J H ))/(^A) 2 , (18) 

and 7rA = 4Tk- Fot models 2 and 3 for the case n = 2 this 
result can be generalized using renormalized perturbation 
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FIG. 6. (Color online) A plot of Xs(Ja)/xa{0) versus 
Jh/4Tk(0) for the models for the data sets shown in Fig. 
□ 



to second order in the renormalized interaction vertices 
for the impurity self-energ y 3 ' 10 ' 11 to give 

/= (2L> 2 +4C/ 1 2 2 + 3Jg)/(7rA) 2 , (19) 

which takes the value I = 10/3 in the regime where Eq. 
(H) holds. 



IV. CONCLUSIONS 

We have shown that we can get very different results 
for the Kondo temperature as a function of the Hund's 
rule coupling Jh depending on how this term is included 
in a model of a magnetic impurity. The physically most 
relevant model for a magnetic impurity with n-fold de- 
generate states in the absence of crystal field or spin-orbit 
splitting should be the one introduced by Yoshimori 6 
which conserves orbital angular momentum about the 



impurity site, or the generalizations of this model con- 
sidered by Mihaly and Zawadowsky 12 , and Yoshimori 
and Zawadowski 7 , where off-diagonal matrix elements 
of the scattering between orbital channels are included, 
subject to the condition of overall conservation of an- 
gular momentum. These matrix elements omitted from 
the Yoshimori model have been shown by Nozieres and 
Blandini 3 - not to affect the results in low temperature 
Fermi liquid regime, nor the derivation of the impurity 
spin and orbital susceptibilities based on the Ward iden- 
tities deduced from conservation of spin and orbital angu- 
lar momentum-. The restriction to models where angular 
momentum is conserved about the impurity site might 
not be appropriate to all situations, such as use of the 
model 3 (given in Eq. to describe the interaction be- 
tween two quantum dots^. We have studied in particular 
the model used by Nevidomskyy and Coleman^ for the 
case n — 2 and, for the smaller values of Jr, verified their 
scaling equation Tj< oc 1/Jh when the Hund's rule cou- 
pling becomes strong enough to lock the individual spins 
to form an effective spin 1 coupled to the conduction elec- 
trons with an effective Kondo coupling J^. We also found 
the point Jh ~ Jk where this scaling crosses over to the 
regime considered by Schrieffer 1 , where J£ = Jk/ti and 
the Tk becomes independent of Jh- 

We find that the low temperature behaviour of all 
the models considered can be described in terms of a 
quasiparticle Hamiltonian with renormalized parameters, 
which, for the case n = 2, can be deduced from an analy- 
sis of the low energy fixed point of the NRG calculation. 
We have shown that, despite the differences between the 
models, in the strong Hund's rule regime, just two inter- 
action parameters are required, a renormalized interac- 
tion U , within each impurity orbital, and a renormalized 
Hund's rule interaction between electrons in different or- 
bitals Jh, which can both be expressed in terms of the 
Kondo temperature, U = 4Tk and Jh = 8Tk/3. 
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